First Estimations of Cosmological Parameters from Boomerang 
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The anisotropy of the cosmic microwave background radiation contains information about the 
contents and history of the universe. We report new limits on cosmological parameters derived from 
the angular power spectrum measured in the first Antarctic flight of the Boomerang experiment. 
Within the framework of inflation-motivated adiabatic cold dark matter models, and using only 
weakly restrictive prior probabilites on the age of the universe and the Hubble expansion parameter 
h, we find that the curvature is consistent with flat and that the primordial fluctuation spectrum 
is consistent with scale invariant, in agreement with the basic inflation paradigm. We find that 
the data prefer a baryon density Q,^h 2 above, though similar to, the estimates from light element 
abundances and big bang nucleosynthesis. When combined with large scale structure observations, 
the Boomerang data provide clear detections of both dark matter and dark energy contributions 
to the total energy density Otot , independent of data from high redshift supernovae. 



The angular power spectrum Ci of temperature 
anisotropy in the cosmic microwave background (CMB) 
is a powerful probe of the content and nature of the uni- 
verse. The DMR instrument on the COBE satellite mea- 
sured Ci for multipoles £ < 20, corresponding to angular 
scales > 7°[Q. Significant experimental effort by many 
groups focusing on smaller angular scales, when com- 
bined J|, ||, Q], led to the Ci estimates in the £ bands 
marked with x's in Figure [IJ which indicate a peak at 
I ~ 200 [gj. It has long been recognized that if Cg can be 
determined with high precision over these angular scales, 
parameters such as the total energy density and baryon 
content of the universe, and the shape of the primordial 
power spectrum of density fluctuations, can be accurately 
measured [[|. The most recently published Boomerang 
angular power spectrum shown in Figure [l] represents 
a qualitative step towards such high precision |7j (here- 
after, B98). 

The data define a strong peak at £ ~ 200. The steep 
drop in power from £ ~ 200 to £ ~ 400 is consistent with 
the structure expected from acoustic oscillations in adi- 
abatic cold dark matter (CDM) models of the universe, 
but is not consistent with the locations and widths of 
peaks expected in the simplest cosmic string, global topo- 



logical defect, and isocurvature perturbation models M. 
The data at higher £ show strong detections which limit 
the height of a second peak, but are consistent with the 
height expected in many CDM models. 

In this paper, we concentrate on determining a set 
of 7 cosmological parameters that characterize a very 
broad class of CDM models by statistically confronting 
the theoretical C^'s with the B98 and DMR data. Sam- 
ple CDM models that fit the data are shown in Figure [|. 
These arc best-fit theoretical models using successively 
more restrictive "prior probabilities" on the parameters. 
A major theme of this paper is to illustrate explicitly 
how inferences that are drawn from the CMB data de- 
pend on the priors that are assumed. Some of these pri- 
ors are quite weakly restrictive and are generally agreed 
upon by most cosmologists, for example that the Uni- 
verse is older than 10 Gyr and that the Hubble constant 
H = 100 h km s _1 Mpc _1 lies between 45 and 90. More 
strongly restrictive priors rely on specific measurements, 
e.g., the HST key project determination of Hq to 10% 
accuracy |9| and the determination of the cosmological 
baryon density, = flbh 2 , to 10% JTofl . 

In j7| , we applied a "medium" set of priors to the B98 
power spectrum to constrain a 6 cosmological parame- 
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FIG. 1: CMB angular power spectra, Ce = £(£ + 
l)(|?£m| 2 )/(27r), where the Ti m are the multipole moments 
of the CMB temperature. The closed (green) circles show 
the B98 data. The magenta crosses are a radical compres- 
sion of all the data prior to B98 into optimal bandpow- 
ers [i]], showing the qualitative improvement provided 

by B98 except in the £ < 20 DMR regime, where the COBE 
data are represented as a single bandpower (open black cir- 
cle). (Note that the B98 and prior CMB points at I = 150 
lie on top of each other.) The smooth curves depict power 
spectra for several maximum likelihood models with different 
priors chosen from Table [E| with (fitot , ut, <-^c, ^a, n s , r c ) as 
follows: PI, short dashed line,(1.3,0. 10,0.80,0. 6,0. 80,0.025); 
P4, dot-dashed line,(l. 15,0.03,0. 17,0.4,0.925,0); P8, short- 
long dashed line,(1.05,0.02,0.06,0. 90,0. 825,0); Pll, solid 
line, (1.0,0. 03,0. 17,0. 70,0.95,0.025). These curves are all rea- 
sonable fits to the B98+COBE data. For comparison, 
we plot a Ho = 68, Qa = 0.7 "concordance model" 
which does not fit (dotted line labelled C), with parameters 
(1.0,0.02,0.12,0.70,1.0,0). 



ter model and found a 95% confidence limit for fl to t of 
0.88 < fitot < 1-12. Row P0 of Table | shows the result 
for our full 7 parameter set with a similar medium prior 
(here taken to be h = 0.65 ±0.1, uj b = 0.019 ±0.006, with 
Gaussian errors for both). As we progress through the 
Table, we show the effect of either weakening or strength- 
ening the prior from this starting point. 

Two of our parameters are fundamental for describ- 
ing the physics of the radiative transport of the CMB 
through the epoch at z ~ 1100, when the photons decou- 
pled from the baryons. These are io b and the CDM den- 
sity uj c = Q c h 2 . The acoustic patterns at decoupling are 
related to the sound-crossing distance at that time, r s , 
which is sensitive to these parameters. We fix the density 
of photons and neutrinos |ll| , which are other important 
constituents at this epoch. The observed B98 patterns 



are also sensitive to the "angular diameter distance" to 
photon decoupling, mapping the z ~ 1100 spatial struc- 
ture to the angular structure, and, through its depen- 
dence on geometry, to fltot, the total energy in units of 
the critical density. When fi t ot < 1 (open models), r s is 
mapped to a small angular scale; when f2 t ot > 1 (closed 
models), r s is mapped to a large angular scale. 

This mapping also depends upon the density asso- 
ciated with a cosmological constant, Q\, and fi m = 
(lu c + ujb)/h 2 , as well as on fi^ = 1 — f2 to t, the den- 
sity associated with curvature. Combinations of Qk/^m 
and Q\/Q m which give the same ratio of angular diam- 
eter distance to sound horizon will give nearly identical 
CMB patterns, resulting in a near degeneracy that is 
broken only at large angular scales where photon trans- 
port through time-varying gravitational potentials plays 
a role. One implication of this is that S1a cannot be well 
determined by our data alone, in spite of the high pre- 
cision of B98. We have paid special attention to such 
near-degeneracies |l2j throughout our analysis. 

The universe reionized sometime between photon de- 
coupling and z ~ 5. This suppresses Ci at small scales 
by a factor e~ 2T °, where r c , our fifth parameter, is the 
optical depth to Thompson scattering from the epoch at 
which the universe reionized to the present. 

Our last two parameters characterize the nature of the 
fluctuations arising in the very early universe, through 
a power law "tilt" n s and an overall amplitude factor 
for the primordial perturbations. The simplest inflation 
models have a nearly scale invariant spectrum character- 
ized by n s 1. Of course, many more variables, and 
even functions, may be needed to specify the primordial 
fluctuations, in particular those describing the possible 
contribution of gravity waves, whose role we have also 
tested For our overall amplitude parameter, we use 
In Cio where Cio is the CMB power in the theoretical spec- 
trum at £ = 10. If we wish to relate the CMB data to 
large scale structure (LSS) observations of the Universe, 
we use In erf as the amplitude parameter, where <j\ is the 
model power in the density fluctuations on the scale of 
clusters of galaxies (8/i _1 Mpc). 

Our adopted parameter space is therefore 
{w b ,w c ,ri t ot,^A,«s,T c ,lnCio}. The amplitude C w 
is a continuous variable, and the rest are discretized 
for the purpose of constructing the model database 
we use to compare data and theory. The number of 
values and coverage are: 15, over 0.1 < f2 to t — 1-5; 
14, over 0.0031 < uj b < 0.2; 10, over 0.03 < lu c < 0.8; 
11, over < Q A < 1.1; 9, over < r c < 0.5; 31, over 
0.5 < n a < 1.5. The spacings in each dimension are 
uneven, designed to concentrate coverage in the regions 
preferred by the data and yet still map the outlying 
regions ju]. Fast CMB transport programs |l^] were 
used to construct our Cg databases. Use was made of 
the angular-diameter distance degeneracy and £-space 
compression to reduce the size and computational 
requirements needed to construct such a database. 

Parameter estimation is an integral part of the B98 
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analysis pipeline, which makes statistically well-defined 
maps and corresponding noise matrices from the time- 
ordered data, from which we compute a set of maxi- 
mum likelihood bandpowers, Cb- The likelihood curva- 
ture matrix Tbb' is calculated to provide error estimates 
including correlations between bandpowers. The curva- 
ture matrix Tbb 1 and the curvature matrix evaluated 
at zero signal, J~ B b'i are use d in the offset-lognormal ap- 
proximation to compute likelihood functions L{x,y) = 
P(Cb\x, y) for each combination of parameters x and y in 
our database. Here x is the value of the parameter we are 
limiting, y specifies the values of the other parameters. 

We multiply the likelihood by our chosen priors, and 
marginalize over the values of the other parameters y, 
including the systematic uncertainties in the beamwidth 
and calibration of the measurement jl6| . This yields the 
marginalized likelihood distribution 

C{x) = P(x\C B ) = J P P rior(x,y)L(x,y)dy. (1) 

For clear detections, central values and la limits for the 
explicit database parameters mentioned above are found 
from the 16%, 50% and 84% integrals of £(x). When 
no clear detection exists, these errors can be misleading, 
so for these cases we shift to likelihood falloffs by e -1 / 2 
from the maximum, or variances about the mean of the 
distribution C(x). The mean and variance are used to 
set the limits on other "auxiliary" parameters such as 
h and fl m , which may be nonlinear combinations of the 
database variables. For good detections the three esti- 
mation methods give very good agreement, and yield 2a 
errors that are roughly twice the la ones generally re- 
ported in this paper. 

We have used this method to estimate parameters, us- 
ing the B98 power spectrum of Figure with the COBE 
bandpowers determined by 0] and a variety of priors. 
The results are summarized in Table |: likelihood func- 
tions for selected parameters and priors are shown in Fig- 
ure ||. 

In the presence of degenerate and ill-constrained com- 
binations of parameters, as with CMB data, the edges 
of the data-base form implicit priors. We have con- 
structed our database such that these effective priors are 
extremely broad. This allows us to probe the dependence 
of our results on individually imposed priors. The choice 
of measure on the space is itself a prior; we have used a 
linear measure in each of our variables jl7| . Sufficiently 
restrictive priors can break parameter degeneracies and 
result in more stringent limits on the cosmological pa- 
rameters. Artificially restrictive databases or priors can 
lead to misleading results; thus, priors should be both 
well motivated and tested for stability. We therefore re- 
gard it as essential that the role of "hidden priors" in any 
choice for Ct database construction be clearly articulated. 

To illustrate the effects of the database structure and 
applied priors, we have plotted likelihood functions found 
using only the database and priors (and no B98 data) in 
Figure 0. These should be compared with those plot- 



ted in Figure |2| which include the B98 constraints, as 
discussed below. We now turn to the results found by 
applying different priors, in the general order of weakest 
to strongest applied priors. 

Our "entire database" analysis prefers closed models 
with very high ujb, as shown in line PI of Table || and in 
Figure |[ The low sound speed of these models couples 
with the closed geometry to fit the peak near I ~ 200. 
These models require very high values of h and u>b, and 
have extremely low ages, so we have mapped out this re- 
gion using a coarse grid. The dual-peaked projected like- 
lihood functions shown are reflections of the the complex- 
ity of the full 9-dimensional likelihood hypersurface. We 
note that parameter combinations that appear to have a 
low probability based on the projected one-dimensional 
limits can fit the data quite well, e.g., the fitot = 1 best- 
fit model of Figure |l|. 

Applying weakly restrictive priors (lines P2-P4 in Ta- 
ble |) moves the result away from the low sound speed 
models, to a regime which is stable upon application of 
more restrictive priors, as shown in panels 1 (top left) 
and 4 (bottom left) of Figure |^. Given their gross conflict 
with many other cosmological tests we do not advocate 
the "entire database" models as representative of the ac- 
tual universe, and we proceed with prior-limited analyses 
below. 

The analysis with weakly restrictive priors (P2-P4) 
finds that the curvature is consistent with flat, while 
favoring slightly closed models. The migration toward 
flk = as more restrictive priors are applied, as shown 
in Table | and in panel 1 (top left) of Figure |^, suggests 
caution in drawing any conclusion from the magnitude of 
the likelihood drop at f2& = 0. In fact, as is evident from 
Figure [l], there are models with f2& = that fit the data 
quite well. The likelihood curve simply indicates that 
there are more models with f2^ < than with £1^ > 
that fit the data well. 

We have taken special care to study the effect on the 
likelihood distributions of choosing a different parame- 
terization of our database. For example, we have inves- 
tigated the likelihoods that result from a finely gridded 
database that uses Q c , fib, and h in place of flk, tOb, 
and uj c to determine flk- This second method, restricted 
to r c = models, uses these different variable choices, 
gridding, and a completely different procedure and code 
which uses maximization of the likelihood over other vari- 
ables rather than marginalization. To compare with this 
second method, we have taken the database used for the 
table and mimicked the effective priors due to the pa- 
rameter limits of the second database. The results found 
using these two parameterizations and codes agree quite 
well - in all cases the likelihood curves shift by at most 
a small fraction of their width. For example, for the ap- 
plied priors of P2 the 95% confidence limits on O tot shift 
from 0.99 < Otot < 1-24 for the method used in the table 
to 0.94 < f2 t ot < 1-27 for the fl c , fib, and h method. Due 
to the very steep slope of the likelihood near flk = 0, 
however, even this small shift changes L(flk = 0)/L max 
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from 0.2 to 0.8. We also find L(n k = 0)/L max w 0.8 
if we use maximization, rather than marginalization, in 
the code used for the table. Additionally, we note that 
a downward shift of about 5% in O to t occurs if the 10 
Gyr age constraint is removed from P2. These points, 
plus the obvious compatibility of the data with the best- 
fit fltot — 1 curve in Fig. |l|, reinforces our conclusion 
that there is no significant evidence in the B98 data for 
non-zero curvature. The only valid conclusion to draw 
from the data that we analyze in this paper is that the 
geometry of the universe is very close to flat. 

The baryon density u>b is also well constrained. While 
our results are higher (~ 3a) than the ujf, estimates from 
light element abundances |l0| , it is most remarkable that 
our entirely independent method yields a result that is so 
close to the BBN value. The scalar spectral index n s is 
very stable once weak priors are applied, and is near the 
value expected from inflation. This weak prior analysis 
does not yield a significant detection of Oa; the Vt c h 2 
results in Table | are suggestive of a detection, but are at 
least in part driven by the weak priors acting on limits 
of the database [|l8| |l^] as shown in Figure [| The values 
of t c are in the range of expectation of the models, but 
there is no clear detection. 

Note that the the weak priors adopt the conservative 
restriction that the age of the universe exceeds 10 Gyr. 
Without this, the weak h prior still allows a contribu- 
tion, albeit reduced, from the high cob, low sound speed, 
low age solution. With the age restriction, this solution 
is eliminated, and the weak BBN prior (u>b < 0.5) does 
not significantly change the constraint: thus, the "weak 
/i+BBN+age" (P4) and "weak /i+age" (P2) rows are es- 
sentially identical. 

In row P4a, we add a "CMB prior" , which is a full 
likelihood analysis of all prior CMB experiments com- 
bined with B98 and DMR, including appropriate fil- 
ter functions, calibration uncertainties, correlations, and 
noise estimates for use in the offset-lognormal approxi- 
mation [^J. As would be expected given the errors we 
compute on the compressed bandpowers of these experi- 
ments in Figure [l] cf. those for B98, this CMB prior only 
slightly modifies the B98-derived parameters, with n s the 
most notable migration. None the less, as much previ- 
ous analysis of the prior heterogeneous CMB datasets 
has shown |Q, reasonably strong cosmological conclu- 
sions could already be drawn on n s and f2tot- Row P4b 
shows results excluding B98, for the weak prior case, 
through our machinery. Though n s and flk have detec- 
tions consistent with the B98 results, no conclusions can 
be drawn on u>b (though the "whole database" analysis 
does pick up the high ojb, Sltot region). We note that if 
t c w is enforced, most variables remain unmoved, but 
n s , which is well-correlated with t c , moves closer to unity: 
for P4, P4a, P4b, we would have n s = 0.97, 1.03, 1.02, 
respectively, and for P5, P5a, P5b, we would have n s — 
0.93, 0.98, 0.98. A prior probability of r c based on ideas of 
early star formation would help to decrease the n s degree 
of freedom. 



The fitot, <-^>b, and n s results are stable to the addition 
of a prior which imposes two constraints derived from 
LSS observations [|lj. The first is an estimate of a\ 
that requires the theory in question to reproduce the lo- 
cal abundance of clusters of galaxies. The second is an 
estimate of a shape parameter for the density power spec- 
trum derived from observations of galaxy clustering p2[ . 
Adding LSS to the weak h and BBN priors (P5, and pan- 
els 2 (top center) and 3 (top right) of Figure ||) breaks 
a degeneracy, yielding a detection of fl\ that is consis- 
tent with "cosmic concordance" models. This also occurs 
when LSS is added to only the prior CMB data (P5b 
and j2l|]). The LSS prior also strengthens the statisti- 
cal significance of the determination of fl c h 2 . Panel 3 
(top right) of Figure ||| shows likelihood contours in the 
Qk = 1 — Qtot vs. plane. Here we have plotted the 
LSS prior (P5), which strongly localizes the contours 
away from the — axis, toward a region that is highly 
consistent with the SNla results [^i). Indeed, treating 
the SNla likelihood as a prior does not change the re- 
sults very much, as indicated by row P12 and P13 of the 
Table, to be compared with rows P5 and Pll, respec- 
tively. 

The use of a strong h prior alone yields results very 
similar to those for the weak h case. The strong BBN 
prior, however, shifts many of the results from the weak 
BBN case. Our data indicate a higher f^/i 2 than BBN, 
and constraining it with the BBN prior shifts the values 
of several parameters, including fl c h 2 , n s , and O m . 
Additional "strong prior" results (P8-P11) are shown in 
Table |, as an exercise in the power of combining other 
constraints with CMB data of this quality. 

A number of the cosmological parameters are highly 
correlated, reflecting weak degeneracies in the broad but 
restricted ^-space range that the B98+DMR data cov- 
ers [[[2] . Some of these degeneracies can be broken with 
data at higher I, as is visually evident in the radically dif- 
ferent behavior of the models of Figure beyond ^ - 600. 
To understand the degeneracies within the context of this 
data, we have explored the structure of the parameter co- 
variance matrix (AyiAt/j), both for the database param- 
eters and the ones derived from them. They add motiva- 
tion for the specific parameter choices we have made [^5| . 
Parameter eigenmodes || [l2| of the covariance matrix, 
found by rotating into principal components, explicitly 
show the combinations of physical database variables 
that give orthogonal error bars. A by-product is a rank- 
ordered set of eigenvalues, which show that for the cur- 
rent B98 data, 3 combinations of the 7 parameters are 
determined to better than 10% p6| . 

We conclude that the B98 data are consistent with the 
predictions of the basic inflationary paradigm: that the 
curvature of the universe is near zero (ilk — 0) and that 
the primordial power spectrum is nearly scale-invariant 
(n s = 1). The slight preference that the current data 
show for closed, rather than open, models is not, we 
believe, a statistically significant indication of non-zero 
curvature. A more conclusive statement awaits further 
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analysis of B98 data, which will increase the precision of 
the measured power spectrum, and/or the results from 
other experiments. 

We measure a strong detection of the baryon density 
0;,/i 2 , a first for determinations of this parameter from 
CMB data. The value that we measure is robust to the 
choice of prior, and is both remarkably close to and signif- 
icantly higher than that given by the observed light ele- 
ment abundances combined with BBN theory. Assuming 
that fitot = 1, we find fl b h 2 = 0.031 ± 0.004. 

Finally, we find that combining the B98 data with our 
relatively weak prior representing LSS observations and 
with our other weak priors on the Hubble constant and 
the age of the universe yields a clear detection of both 
non-baryonic matter {£l c h? — O.OITIqqq^) and dark en- 
ergy (Q\ = 0.661q'q9) contributions to the total energy 
density in the universe. The amount of dark energy that 
we measure is robust to the inclusion of a prior on f2 to t 
(shifting to Q\ = 0.67^o;gg for Q to t — 1), and to the 
inclusion of the prior likelihood given by observations 
of high-redshift SNla (shifting to fi A = 0.69±g;g| when 
both the fitot = 1 and the SNla priors are included). The 



perfect concordance between the completely independent 
detections of J7a from the CMB+LSS data and from the 
SNla data is powerful support for the notion that the 
universe is currently dominated by precisely the amount 
of dark energy necessary to provide zero curvature. 

The analysis presented here and in Q makes use of 
only a small fraction of the data obtained during the first 
Antarctic flight of Boomerang. Work now in progress 
will increase the precision of the power spectrum from 
I = 50 to 600, and extend the measurements to smaller 
angular scales, allowing yet more precise determinations 
of several of the cosmological parameters. 
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They have little effect over the range of £'s that B98 is 
most sensitive to, but could have an important impact 
on the amplitude relative to COBE. To test the role that 
GW induced anisotropics would play, we have adopted 
the model used by j21| : for n s < 1, we set n t — n s — 1 
and for n s > 1, we allow no GW contribution. This 
presents a fixed alternative, reasonably motivated by in- 
flation, without introducing new parameters. We have 
found that there is a negligible effect on the parameter 
determinations in Table [ft there is only a very slight mi- 
gration upward in n s . 

The specific values of the database variables used for this 
analysis are: (lu c = 0.03, 0.06, 0.12, 0.17, 0.22, 0.27, 0.33, 
0.40, 0.55, 0.8), (w h = 0.003125, 0.00625, 0.0125, 0.0175, 
0.020, 0.025, 0.030, 0.035, 0.04, 0.05, 0.075, 0.10, 0.15, 
0.2), (Ov =0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 
1.1), (Q k =0.9, 0.7, 0.5, 0.3, 0.2, 0.15, 0.1, 0.05, 0, -0.05, 
-0.1, -0.15, -0.2, -0.3, -0.5), (n s =1.5, 1.45, 1.4, 1.35, 1.3, 
1.25, 1.2, 1.175, 1.15, 1.125, 1.1, 1.075, 1.05, 1.025, 1.0, 
0.975, 0.95, 0.925, 0.9, 0.875, 0.85, 0.825, 0.8, 0.775, 0.75, 
0.725, 0.7, 0.65, 0.6, 0.55, 0.5), (r c =0, 0.025, 0.05, 0.075, 
0.1, 0.15, 0.2, 0.3, 0.5) 

U. Seljak and M. Zaldarriaga, Astrophys. J. 46 9, 



437 (1996L http://www.sns.ias.edu/~matiasz/CMB- 



O'Meara, N. Suzuki and D. Lubin, submitted to Physica 



Scripta, istro- pli/0001318| (2000). 
The Ce spectra with massive neutrinos are quite similar 
to those without, and current data, including B98, will 
not be able to strongly constrain the value. When the 
LSS prior is added to the CMB data, however, the com- 
bination is quite powerful, e.g., |2l| ], 

G. Efstathiou and J. R. Bond, Mon. Not. R. Astron. Soc. 
304, 75 (1999). 

Gravity waves (GW) can induce CMB anisotropy, and 
could have a separate tilt, n t , and an overall amplitude. 



FAST / cmbfast .litiii 



A. Lewis, A. Challinor and A. 



Lasenby, |astro-ph/9911177| (1999). 
[16] Apart from the 7 stated database parameters, we have 
allowed for an estimated 10% uncertainty in the calibra- 
tion and the beam, which we determine simultaneously 
with the overall amplitude Cio, by relaxing to the maxi- 
mum likelihood value in these variables. We then deter- 
mine the Fisher error matrix, assume that the variables 
are log-normally distributed, and evaluate a correction to 
the likelihoods appropriate for marginalization over these 
"intrinsic" variables. Including the marginalization cor- 
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TABLE I: Results of parameter extraction using successively more restrictive priors. The confidence intervals are la, evaluated 
using methods described in the text. The 2a errors are approximately double the la values quoted in most cases; upper limits 
are quoted at 2a. The quoted values are reported after marginalizing over all other parameters. Note that these combinations 
are not, and should not be, the parameters of the "maximum likelihood" best-fit models of Figure jjj. The weak h and BBN 
(fii,/) 2 ) priors are tophat functions (uniform priors) and both include an additional age > 10 Gyr prior. The strong priors are 
Gaussians with the stated la error, and also have weak constraints imposed on the other variables. P0 is the medium h + 
BBN prior used in |Q and described in the text. The LSS priors are combinations of Gaussians and tophats The SN1 

prior (P12, P13) includes LSS as well the SNla likelihood shown in panel 3 of Fig. |ij P4a and P5a show the small effect of 
including prior CMB data in our B98+DMR analysis; these should be contrasted with P4b and P5b, the case of prior CMB 
data alone. Columns 1-5 (fitot to Oa) are predominantly driven by the CMB data, except for fi^/i 2 and fif, when the strong 
BBN prior (P7-P9) is applied. Most of the values in columns 6-10 (r c to Age) are influenced by the structure of the parameter 
space and should not be interpreted as CMB-driven constraints; exceptions are the fi c /i 2 and fiA results when the LSS prior is 
applied. An equivalent table that includes an inflation-inspired gravity wave induced contribution to the anisotropy [I13J yields 
remarkably similar parameters and errors. 

Priors fitot Q,bh 2 n s fif, fiA t c fi c /i 2 fi m h Age 

po: Medium /i+bbn 1.078:81 o.o3og;88t i.oo8:8l 0.088:81 o.378:i 0.128:^ o.258:J8 0.728:! o.63g : 8l n.9\i 



pi 



Whole Database OlbVie 0.1008:8t5 0.888:51 0.108:81 0.538:1? 0.228:11 ... o.8i8:ll i.088:i 7.82.I 



P2 
P3 
P4 



Weak h (0.45 <h< 0.90)+age 1.158:58 0.0368:881 1-048:58 0.1lg$ < 0.83 0.218:11 0.24g:g| 0.848:1 0.58g:18 12-7? 1 



2.1 



Weak BBN {fl b h 2 < 0.05)+age 1.168:18 0.0358:881 l-038:io 0.168:81 < 0.83 0.218:11 0.198:o9 0.928:i 0.528:11 14.6l:| 

Weak fe+BBN+age i.i58:58 0.0368:881 i-048:58 o.n8:81 < o.83 0.218:11 o.248:81 o.848:I 0.588:18 12.7II 



P4a: Weak and prior CMB 1-018:81 0.03l8:88e 1-068:58 0.108:8! < 0.79 0.248:1? 0.188:8e 0.648:i 0.598:11 13.41:8 

P4b NO B98: Weak and prior CMB LOggjg 0.0248:81s 1-148:11 0-088:88 < 0-80 0.298:18 0-2lg;gl 1 7lgJ| 0.608:11 12.9|:g 

P5: LSS & Weak /i+BBN+age 1.128:8? 0.0348:881 0.998:5s 0-108:81 0.668:89 0.19g:?i 0.148:81 0.488:11 0.608:11 14.51:8 

P5a: lss & Weak and prior cmb i.028:81 0.0308:880 1-058:58 0.098:8! 0.478:51 0.228:11 0.168:8! 0.578:18 0.608:11 13.81:? 

P5b NO B98: LSS & Weak and CMB 1.008:8? 0-0288:811 1-088:11 0.088:88 0.588:1? 0.268:11 0.148:81 0.448:11 0.638:11 13.81:? 

,0.07 n flQfiO.005 1 ncO.09 n aqO.03 ,noo n on0.19 n ocO.OS n -71 0.27 n ccO.07 11 fil.4 



P6 
P7 
P8 
P9 



strong h (h = 0.71 ± 0.O8) i-098:8e 0.0368:881 i-058:81 0.088:83 < 0.82 0.208:1! o.26g:?§ 0.718:1? 0.668:8? 11.6J 

strong bbn (fi 6 h 2 =o.oi9 ± 0.002) i.iogiSS o.o2i8:88l o.858:8? 0.078:81 o.798:88 0.098:5? o.os8:8I o.38g:li 0-548:18 i7.7l:l 

strong /i+bbn i.048:81 o.o2ig : 881 o.878:8? 0.058:81 o.758:M o.os8:J| o.o98:8l 0-288:11 o.6s8:88 i5.2l : | 

LSS fc Strong fe+BBN 1.048:81 0-0228:881 0-928:88 0-058:81 0.66g : g? 0-088: | OTfgjf 0.398:8? 0-64g : 81 14-Ojj 

pio: n tot = 1 & weak /i+age i o.o3i8:88! 0.998:8? 0.068:81 < o.78 o.io8:5? 0.278:87 o.578:ll o.748:81 io.98:| 

pii: fitot = i & lss & Weak i o.o30g:88t 0-968:8e 0-058:81 o.678:8o o.o98:J? o.i88:8i 0.328:8! 0.798:8! n.78:! 



P12 

P13 



lss & Weak & sNia i.088:8! 0.0348:88! 1.028:81 o.o88:8i 0.728:8! 0.238:!? o.is8:81 0.378:8? 0.708:88 i3.3l:| 

fitot = 1 & lss & weak & SNia i o.o3og;ggt o.97g:g£ 0.058:81 o.698:8! o.io8:5? 0.188:8? 0.318:8! o.8i8:81 n-68:l 



rection makes little difference. We have also marginalized 
over bins that were used in creating the power spectrum 
but not in the analysis, since they are correlated. 
[17] The choice of measure is not important for strong local- 
ized peaks, but can potentially affect limits on poorly 
constrained variables and on those with complex likeli- 
hood functions. One can argue for logarithmic measures 
in C10 (as we have used) and in ujb and u) c (which we 
have not used), and there are certainly philosophical al- 
ternatives to linear measures in fitot and fiA- Consider 
what happens when we turn the "whole database" w;, 
likelihood curve of Figure ^ into a probability curve if 
we adopt a logarithmic rather than linear measure: the 
anomalous peak at 0.1 drops below the "cosmological 
peak" at 0.03; once weak priors are adopted, the 0.03 
peak is all that is left and it is very stable to changing 
the measure. Changing measures usually moves peaks a 
small fraction of a a, although the amount does depend 
upon relationships to correlated variables with large er- 
rors. The discreteness of our database is also a restric- 
tion on how accurately we can localize peaks. For exam- 
ple, our finest gridding in fitot is 0.05 from 0.8 to 1.2, 



hence accurate localization better than half this spacing 
should not be expected. When projections are made, the 
available volume of models leads to effective priors as 
well@|l9|. 

[18] The weak prior by itself actually focuses u) c about 0.22, 
dropping to either side because of h and age restrictions. 
Our data do constrain lo c further, but not enough to claim 
a CMB determination beyond the prior until the LSS 
prior is included. 

[19] fiA and fitot have a prior probability dropping as fitot 
drops and tt\ rises just because fi m is positive. There is 
a physical effect that also favors the closed models when 
CMB is added. As u>b varies, the sound speed lowers, the 
peak moves to higher t, but can be mapped back to our 
observed £ p k by judiciously choosing an fitot > 1. fiA > 
moves the peak to lower I which fitot < 1 can also move 
back to the observed £ p k, but it is a smaller effect. If we 
had allowed fiA < 0, closed models could have done the 
same, further favoring fitot > 1 because of the volume of 
models available. 

[20] e.g., M. White et ai, Mon. Not. R. Astron. Soc. 283, 
ppl07, (1996); K. Ganga et al, Astrophys. J. 484, 7 



7 



(1997); C. Lineweaver, Astrophys. J. 505, 69 (1998); ref. 
Q; G. Efsthathiou et al, Mon. Not. R. Astron. Soc. 
303, pp47, (1998); S. Dodelson and L. Knox, Phys. Rev. 
Lett. 84, 3523 (2000); A. Melchiorri et al, Ap. J. Let- 
ters, 536, L63, (2000); M. Tegmark and M . Zaldarriaga, 
Astrophys. J., in press, |a.stro-ph/000209l| (2000); M. Le 



Dour et al. submitted to Astron. and Astrophys., astro 



ph/0004282 (2000) 



[21] J.R. Bond and A.H. Jaffe. Phil. T rans. R. Soc. London 
357, 57 (1999), |astro-ph/9809043| . 

[22] The LSS prior is a slight modification of the one used 
in j pi| . <78f257i 56 =0.55l'Q2'^'o8 ^ s assumed to be dis- 
tributed as a Gaussian smeared by a tophat distribu- 
tion, with the first error indicating the 1-a error on 
the Gaussian, and the second indicating the extent of 
the tophat about the mean. The constraint from power 
spectrum shapes involves a combination of spectrum 
tilt, n s — 1, and a "scaling shape parameter" T ~ 
O^e-^d+^W^-O-oe) which is related to the 
horizon scale when the Universe passed from relativis- 
tic to matter dominance: F + (n s — 1)/2=0.22^'q4'^'qy. 
Both priors were designed to generously encompass the 
observations, and so are "weak" to "medium" rather than 
"strong" , in the sense of Table |l[ 

[23] The contours plotted at C/C m ax = 
exp[— {2.30, 6.17, 11. 8}/2] provide rough indicators 



of 1, 2, and 3cr |J. 

[24] Riess et al, Astron. J. 116, 1009 (1998); S. Perlmutter 
et al, Astrophys. J. 517, 565 (1999). 

[25] Here are some sample correlation coefficients for the weak 
fr±BBN case of Table [j| it is relatively small between u>i, 
and other database variables but between Qb and h it 
is 86%. Similarly, as is evident from the contour map in 
Figure [], flk and D.a are correlated only at the 41% level, 
whereas Q m and Qa are correlated at the 96% level. Thus, 
for CMB work it is advantageous to use Qk as a variable 
rather than Q m , and hence this is what we plotted in 
Figure |^ rather than the more recognizable Q m -QA plot. 
lo c and Qa have a 75% correlation, not surprising in view 
of that for Q m . 

[26] For the P4 case, the best determined eigenmode (to 
±0.03) is a combination of slope, amplitude and Thomp- 
son depth; next (to ±0.04) is predominantly Qk, with a 
judicious negative contribution from Qa, a combination 
orthogonal to the angular diameter distance degeneracy; 
the third eigenmode (to ±0.09) is mostly ujb, with a lit- 
tle contribution from all other variables. The next three 
combinations are determined to ±0.14. the worst (±0.4) 
combination is one of uj c and Qa- Similar coefficients and 
accuracies hold for other priors, except for distortions in 
the strong BBN prior case. 
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FIG. 2: Likelihood functions for a subset of the priors used in Table |l[ Panel 1 (top left) shows the likelihood for = 1 — fitot ; 
the full-database (PI, dotted line) prefers closed models, but reasonable priors (P2, dashed blue line; P4, solid blue line; PO, 
dot-dash red line; note that P2 and P4 lie on top of one another in every panel in this plot but are distinct in Figure ^) 
progressively move toward = 0. We caution the reader against agressively interpreting any 2<r effects. Likelihood curves for 
Qa are shown in panel 2 (top center). In panels 2 and 4-6, the cases and line types are as in panel 1, except that dot-dashed 
now denotes the weak+LSS prior, P5. With weak priors applied, there is no significant detection of S7a (P2 and P4, overlapping 
as solid blue line in all remaining C(x) panels). Only by adding the LSS prior is localized away from zero (P5, red dot-dash 
in all remaining C(x) panels). Panel 3 (top right) shows the contour plot of Qk and Q,a, for which the first two panels are 
projections to one axis. The bold diagonal black lines mark fl m =l and f2 m =0. The blue contours are those found with the 
weak prior (P4), plotted at 1, 2, and 3ap3|. Red contours are similarly plotted for the weak+LSS prior (P5). SNla constraints 
are plotted as the lighter (black) smooth contours, and are consistent with the CMB contours at the la level. Panel 4 (bottom 
left) shows the contours for u>b\ the full database analysis results in a bimodal distribution with the higher peak concentrated 
at very high values. These high ujb models are eliminated by the application of a weak h prior or weak BBN prior (P2 and 
P4, overlapping as blue here). Panel 5 (bottom center) shows a localization of u) c for the weak h and BBN prior cases, but 
this is partially due to the effect of the database structure coupling to the h and age priors. Only the LSS prior (P5, red 
dot-dash) allows the CMB to significantly constrain ui c . Panel 6 (bottom right) shows good localization and consistency in 
the n s determination once any priors are applied. The inflation-motivated fi t ot = l priors (P10, Pll) give very similar curves 
localized around unity. See Figure )a to see the effects of the database and priors on these curves. 
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FIG. 3: Likelihood functions similar to those in Figure but computed without using the B98 data. These curves show the 
effect of the database constraints and applied priors alone. The identification of the curves is the same as in Figure j| with the 
addition of the dotted magenta curve in panels 2-6, which shows the likelihood given weak priors and the COBE DMR data. In 
panel 3, only the la (red) contour is shown for the prior only and prior+DMR cases, while 1, 2 and 3<r (light black) contours 
are shown for SNla. The curves for P2 (solid blue) and P4 (dashed blue) are slightly separated in this figure, in contrast to 
Figure ^, where they overlapped. Of particular interest here are the slope induced across flk, the slight localization of Q c h 2 
with the weak priors, and the significant localization of Q c h 2 and n s with just weak+DMR+LSS (dotted magenta). 



